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ABSTRACT 

. Context. Amplitudes and phases of the light variation of a pulsating star in various photometric passbands contain information about 

' geometry of observed modes. Because oscillation spectra of early B-type main sequence stars do not exhibit regular patterns, these 

observables are very often the only ones from which mode identification can be derived. Moreover, these data can yield valuable 
constraints on mean stellar parameters, subphotospheric convection, microphysics and atmospheres. 

Aims. We study all possible sources of inaccuracy in theoretical values of the photometric observables, i.e. amplitude ratios and phase 
differences, of early B-type main sequence pulsators. Here, we discuss effects of parameters coming from both models of stellar 
' atmospheres and linear nonadiabatic theory of stellar pulsation. In particular, we evaluate for the first time the effect of the departure 

from the LTE approximation. 

Methods. The photometric amplitudes and phases are calculated from the semi-analytical formula for the light variation. The atmo- 
spheric input comes from line-blanketed, LTE and NLTE plane-parallel, hydrostatic models. To compute the limb-darkening coef- 
' ficients for NLTE models, we use the Least-Square Method taking into account the accuracy of the flux conservation. The linear 

. , nonadiabatic stellar pulsations were computed by means of Dziembowski code. We consider the OPAL and OP opacity tables and 

. two determinations of the solar chemical mixture: GN93 and AGSS09. 

{ZLh' Results. We present effects of NLTE atmospheres, chemical composition and opacities on theoretical values of the photometric ob- 

servables of early B-type pulsators. To this end, we compute tables with the passband fluxes, flux derivatives over effective temperature 
and gravity as well as the non-linear limb-darkening coefficients in 12 most often used passbands, i.e. in the Stromgem system, uvby, 
and in the Johnson-Cousins-Glass system, UBVRIJHK. We make these tables public available at the Wroclaw HELAS Web page *. 



o 

^ ' Conclusions. In the case of radial modes, effects of uncertainties in model atmospheres are much smaller than those coming from 

' the pulsation theory. In turn, for nonradial modes effects of NLTE become more important and they are most significant for dipole 

modes. Therefore, if one wants to construct an accurate seismic model and gain more constraints, all inaccuracies in the photometric 
observables have to be allowed for. 
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O ' "I ■ Introduction Stobie (WJ9) and Stamford & Watson (1981') expanded it for the 

light variation in the photometric passbands. They showed that 

a> Data on the photometric and spectroscopic vanations of a pul- ^^^^^ ^^^^ ^iff^^g^j ^^j^^^ ^ j^^^j^^ i„ separated parts on 

sating star bring information on frequencies of excited modes ^j^^ amplitude ratio vs. phase difference diagrams. Subsequently, 

and their geometry The latter property is particularly impor- ^^^^ ^^^^^^ ^^^^ ^ppji^^ ^^^^^^ ^yp^^ pulsating stars 

tant If oscillation spectra are sparse and lack equidistant pat- ^^^^^^ , ^988 , ^ugier, Dziembowski & Pamyatnykh ( fT994b 

terns. This is the case for main sequence pulsators of early B ij^proved the method by including nonadiabatic effects in calcu- 

spectral type, i.e. p Cephei vai'iables. Additionally, using data j^jj^^^ ^ ^ephei stars. Effects of rotation on photomet- 

^ on the photometi-ic amplitudes and phases one can consti-ain observables were studied by Daszyriska-Daszkiewicz et al. 

mean stellar parameters, test efficiency of convection transport ^j^^^ frequency modes and by Townsend (|2003j and 
■ ■ ■ in the case of cooler pulsators, like 6 Scuti stars, (Daszyriska- Daszyriska-Daszkiewicz, Dziembowski & Pamyatnykh ("SO^ 
Daszkiewicz, Dziembowski, Pamyatnykh 2003) or test opacity long-period g-modes. To compute values of the photomet- 
data in the case of B-type pulsators (Daszyriska-Daszkiewicz, ^-^ amplitudes and phases we need two kinds of input. The first 
Dziembowski, Pamyatnykh 2005). comes from stellar model atmospheres and the second one 
Ones of the most popular tools to identify a pulsation mode f^.^^ computations of stellar pulsation. Both data contain various 
are the amplitude ratios and phase differences in various photo- sources of uncertainties 
metric passbands. In the zero-rotation approximation, these ob- 
servables depend on the mode degree, f, but are independent of The goal of this paper is to examine for the first time an in- 
the azimuthal order, m, the intrinsic mode amplitude, s, and the fluence of NLTE effects on theoretical values of the photometric 
inclination angle, /. amplitude ratios and phase differences for early B-type main se- 
The semi-analytical expression for the bolometric light vari- quence pulsators. We check also effects of metallicity and micro- 
ation was formulated by Dziembowski (119771) . Then, Balona & turbulent velocity in the atmosphere. Subsequently, we compare 
these effects with uncertainties coming from the linear nonadi- 



http://helas.astro.uni.wroc.pl abatic theory of stellar pulsation. As an example, we consider 
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main sequence models with a mass of 10 Mq and low degree 
modes with £-0, 1, 2. Here, we neglect all effects of rotation on 
pulsation. 

The structure of the paper is the following. In Section 2, 
we recall the linear formula for the pulsation complex ampli- 
tude in a photometric band. Sections contains description of 
NLTE model atmospheres and results of our computations of the 
band fluxes, corresponding flux derivatives over effective tem- 
perature and gravity, and the nonlinear limb darkening coeffi- 
cients, in 12 photometric passbands: uvbyUBVRIJHK. Tables 
with these data can be downloaded from the Wroclaw HELAS 
web page. Moreover, we study effects of temperature, gravity, 
NLTE, atmospheric metallicity and microturbulent velocity on 
the above mentioned quantities. How these atmospheric uncer- 
tainties translate into the pulsation photometric observables, i.e. 
amplitude ratios and phase differences, is discussed in Section 4. 
Inaccuracies in the photometric observables connected with the 
linear nonadiabatic theory of stellar pulsation are presented in 
Sections. We end with conclusions in Section 6. 

2. Light variation due to stellar pulsation 

Stellar pulsations cause changes of temperature, normal to the 
surface element and pressure. If pulsations are linear and all ef- 
fects of rotation on pulsation can be ignored, then the total am- 
plitude of the light variation in the passband x can be written 
in the following complex form (Daszynska-Daszkiewicz et al. 
I2002] l: 



(1) 



where e is the intrinsic mode amplitude, Y"' - the spherical har- 
monic and / - the inclination angle. The amplitude is given by 
abs(Jl^) and phase by arg(Jl^). The D^^^-f product stands for 
temperature changes, where 



I tjlogCT^I) 
4 d log Teff 



(2) 



;Ff is the flux in the passband x and / is the nonadiabatic com- 
plex parameter describing the amplitude of the radiative flux per- 
turbation to the radial displacement at the photosphere level 



Re{sfY"i'{0,ifi)e- 



Geometrical term, D2/, is given by 

L»2/ = (2 + ^)(l-^), 

and the pressure term, Di by 



D^^f = - 2 + 



GM j dlogg 



M is the disc averaging factor defined by 



Jo 



(3) 



(4) 



(5) 



(6) 



where h^fj) is the limb darkening law, Pf is the Legendre poly- 
nomial and ju is a cosine of the angle between a line of sight and 
the emergent intensity. Remaining parameters have their usual 
meaning. 

In the above expressions, we can distinguish two sorts of in- 
put parameters needed to compute theoretical values of the pho- 
tometric amplitudes and phases. The first input is derived from 



models of stellar atmospheres and these are the flux derivatives 
over effective temperature and gravity (Eq.|2] and |5]), as well as 
limb-darkening and its derivatives (Eq.|2l|5]and|6]). The second 
input comes from the nonadiabatic theory of stellar pulsation and 
this is the /-parameter (Eq.[T]and[3]). 

In this paper, we used the Warsaw-New Jersey 
evolutionary code and the linear nonadiabatic pul- 
sation code of Dziembowski (1977). We considered 
opacity tables from OPAL (Iglesias & Rogers 1996) and 
OP (Seaton 2005) projects, and two determinations of the solar 
chemical composition: GN93 by Grevesse & Noels 1993 and 
AGSS09, a recent one by Asplund et al. 2009. As for models 
of stellar atmospheres, we considered Kurucz models (Kurucz 
120041) and TLUSTY models (Lanz & Hubenv l2007l l. 

3. NLTE line-blanketed model atmospheres 

The most widely used models of stellar atmospheres are line 
blanketed, plane-parallel, hydrostatic models of Kurucz (2004) 
computed within an approximation of local thermodynamic 
equilibrium (LTE). However, in the case of atmospheres of early 
B-type stars, effects of the departure from LTE and a proper 
treatment of line opacity become important. Escalation of the 
quality of spectro-photometric observations calls for a need of 
high resolution and accurate models of stellar atmospheres. 

A grid of non-LTE (NLTE) model atmospheres were 
computed by Lanz & Hubeny (120031 ) and Lanz & Hubeny 
(2007). These are metal line-blanketed, plane-parallel, hydro- 
static model atmospheres of O-type stars (OSTAR2002) and of 
early B-type stars (BSTAR2006), respectively. In their compu- 
tations, they adopted a solar chemical mixture by Grevesse & 
Sauval (119981) . helium to hydrogen abundance of He/H=0.1 by 
number and two values of the microturbulent velocity, ^,=2 and 
10 km/s. Lanz & Hubeny (2003 2007) showed that in the case 
of OB stars, a neglect of NLTE effects causes differences not 
only in spectral lines but also in the continuum level. In the 
near ultraviolet (the Balmer continuum), the LTE fluxes are up 
to 10% higher than the NLTE counterparts. In turn, in the far and 
extreme ultraviolet (the Lyman continuum), the LTE fluxes are 
lower than the NLTE ones. For more details see Lanz & Hubeny 
(I2003II2007I ). 

The grid of the OSTAR2002 models were computed for 12 
values of effective temperatures appropriate to O-type stars, i.e. 
between 27500 and 55000 K with a step of 2500 K, 8 values of 
surface gravities from logg=3.0 to 4.75 with a step of 0.25 dex 
and one microturbulent velocity, ^/=10 km/s. Moreover, 10 val- 
ues of the atmospheric metallicity were considered, (Z/Zoj^t^-l, 
1, 0.5, 0.2, 0.1, 1/30, 1/50, 1/100, 1/1000, 0.0, where Z is the 
metal abundance by mass and Zq is the solar value. 

The grid of the BSTAR2006 models contains 16 values of 
effective temperatures between 15000 and 30000 K and a step of 
1000 K, 13 values of surface gravities in the range from log^ = 
1.75 to 4.75 dex and a 0.25 dex step. In this case, 6 values of 
metallicity were considered, {Z/Zq).^u^=2, 1, 0.5, 0.2, 0.1, 0.0, 
and two values of the microturbulent velocities, f,=2 km/s and 
^,=10 km/ s . Models with ^, = 1 km/s were computed for two sets 
of chemical mixture but only for B-type supergiants (logg < 
3.0). The first mixture was the same as for ^,=2 km/s, i.e. the 
solar composition, and the second one was enriched in helium 
and nitrogen, and depleted in carbon. 

The lower limit for gravity at a given effective tempera- 
ture was determined approximately by the Eddington limit. The 
TLUSTY code used by Hubeny and Lanz (1995) becomes un- 
stable near this limit. 
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3.1. The passband fluxes and their derivatives 

We computed the fluxes, Tx, for the BSTAR2006 NLTE mod- 
els in 12 commonly used photometric passbands, i.e. in the 
Stromgren (uvby) and Johnson-Cousins-Glass (UBVRIJHK) 
systems, according to the formula 



/ T{A)S{A)dA 



(7) 



where S (A) is the response function of the passband x, adopted 
from the Asiago Database on Photometric Systems (Moro & 
Munari 120001) . The integral in Eq. (|7]l is computed in the wave- 
length range from A[ to A2 where S (A) has non-zero values. 

We considered the whole grid of effective temperature, Teff , 
gravity, log g, metallicity, (Z/ZQ)atm, and microturbulent veloc- 
ity, for the BSTAR2006 models. Tab le[T] summarizes our re- 
sults. The file names are coded by values of (Z/ZQ)atm and We 
have adopted the same designations as Lanz &Hubeny ( 120071 ). 
i.e. BC, BG, BL, BS, BT, BZ denote {Z/Zq)^i^=2, 1, 0.5, 0.2, 
0. 1 and 0, respectively, whereas v2 and vlO correspond to ^, = 2 
km/s and 10 km/s, respectively. The index CN marks a model 
enriched in helium and nitrogen, and depleted in carbon. Each 
file contains the following columns: line number, effective tem- 
perature, Teff, logarithm of the surface gravity, log g, metallicity, 
(Z/Z0)atm, microturbulent velocity, and the logarithmic flux, 
log;Fi, in the uvbyUBVRIJHK passhands. 

The most important atmospheric parameters in the expres- 
sion for the brightness variation of a pulsating star are the flux 
derivatives over effective temperature, T^s, and gravity, log g. In 
a given photometric pasband x, they are defined as 



d log Teff 



and 



0! 'a - 



dlogg 



Values of these flux derivatives depend not only on effective tem- 
perature and gravity but also on the metallicity, microturbulent 
velocity and the departure from LTE in the star's atmosphere. 

Similarly to the passband fluxes, the derivatives were com- 
puted in the whole range of parameters of the BSTAR2006 
models. Names of the derivative tables are given in Table[T] 
Each file contains the following columns: line number, effec- 
tive temperature, Teff, logarithm of the surface gravity, log^, 
metallicity, {Z/ZQ)ntm, microturbulent velocity, and the flux 
derivatives over effective temperature, aj, and gravity, a^, in the 
uvbyUBVRIJHK passhands. 

In Fig. 1, we show the NLTE flux derivatives as a function of 
temperature for two Stromgren passbands uy and three values of 
gravity, log g = 3.5, 4.0, 4.5. In the left panel, we plot the temper- 
ature derivative, ax, and in the right one the gravity derivative, 
ag. We assumed the solar metallicity, (Z/Zo)atm=l and the mi- 
croturbulent velocity of ^,-2 km/s. As we can see the lower the 
gravity the larger values of aj and the lower values of ag. 

The wavelength dependence of the flux derivatives is pre- 
sented in Fig. 2. We considered the central wavelengths, Ac, of 
the UBVRIJHK passbands and a model with T^ff = 20000 K 
and log^ =4.0. The most steep derivatives are in the ultraviolet 
filter and then the absolute values of aj and Ug decrease quite 
rapidly with A^. This is because the early B-type stars emit the 
most amount of their energy in the UV wavelength range. In 
Fig. 2, we compare also the NLTE derivatives with the LTE ones. 
The largest difference can be seen again in the U passband. 



Effect of NLTE on the flux derivatives as a function of ef- 
fective temperature for the Stromgren uy passbands is presented 
in Fig. 3. In general, differences between the LTE and NLTE 
derivatives are relatively small and they increase with effective 
temperature, especially in the case of ag in the u passband. In the 
case of both model atmospheres, the values of aj in the u pass- 
band change only slightly with the temperature, whereas in the y 
passband they are larger for higher T^ff. The gravity derivatives 
get more steep with the higher temperature, except ag(Tss) for 
NLTE models, which is a decreasing function up to Teff 24000 
K and then it increases. 



3.2. Limb darkening 

Knowledge of a distribution of the specific intensity over the 
stellar disc is crucial in many fields of astrophysics. In the case of 
pulsating stars, Umb-darkening function and its derivatives over 
effective temperature and gravity, are the second atmospheric 
data needed to compute theoretical values of the photometric 
amplitudes and phases. 

A usage of analytical formulae for the specific intensity dis- 
tribution saves enormously the computation time. The first limb- 
darkening law was proposed by Milne ( 1921 ) in a linear form. 
With the development of models of stellar atmospheres, it turned 
out that this approximation has a poor accuracy. Thereafter, 
Klinglesmith & Sobieski ( 1970) published a logarithmic law for 
early type stars. A quadratic formula for the limb-darkening law 
was proposed by Manduca at al. ( 1977 ) and Wade & Rucinski 
(1985). For stars hotter than 8500 K, Diaz-Cordoves & Gimenez 
( 119921 ) suggested a square root law. 

However, all these above formulas were not adequate for 
all types of stars. A more general law was proposed by Claret 
(120001) in the following non-linear form 



m, 1) 



(8) 



k=\ 



where I{A,iS) is a specific intensity at the wavelength /i, yU is a 
cosine of the angle between a line of sight and the emergent 
intensity, and I{A, 1) is the value at the center of the stellar disc. 
Limb darkening coefficients, are determined to reproduce the 
model intensity distribution and to conserve the flux with a high 
accuracy in the whole range of effective temperatures and gravi- 
ties. 

These coefficients were calculated for LTE model atmo- 
spheres in many photometric systems and a wide range of ef- 
fective temperatures, gravities, metallicities and for several val- 
ues of the microturbulent velocities (Claret 2000, 2003, Claret 
& Hauschildt 2003, Claret 2008). A comparison of different nu- 
merical methods of computations of limb darkening coefficients 
(LDC) were discussed in detail by Diaz-Cordoves et al. (119951) . 
Claret (l2000l l. Heyrovsky (2007) and Claret (2008ll. 

In this paper, we determine non-linear LDC for metal line- 
blanketed, NLTE, plane-parallel, hydrostatic model atmospheres 
of Lanz & Hubeny (2007). In the first step, we computed 
monochromatic specific intensities for 20 equally separated 
points of in the range of (0.001 , 1 ) and in the wavelength range 
of (2950,26500) A, for all parameters of the BSTAR2006 atmo- 
spheres. All calculations were performed using the SYNSPEC 
program (Hubeny at al. 1985 Hubeny & Lanz 2000). This code 
is intended to compute specific intensities and fluxes from the 
model atmosphere input, which we took from NLTE models de- 
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Table 1. Tables of the passband fluxes, their derivatives over effective temperature and gravity, and non-linear limb darkening 
coefficients, for the BSTAR2006 NLTE models in the uvbyU BV RIJ H K passhands. 



fluxes 


derivatives 


limb darkening 
coefficients 


range of T^g [K] 


range of log g [dex] 




^, [km/s] 


mixture 


flux_BCv2 


der.BCv2 


LDC_BCv2 


15000-30000 


1.75-4.75 


2.0 


2 


GS98 


flux_BCvlO 


der.BCvlO 


LDC_BCvlO 


15000-30000 


1.75-3.00 


2.0 


10 


GS98 


ffux_BCvfOCN 


der.BCvlOCN 


LDC_BCvlOCN 


15000-30000 


1.75-3.00 


2.0 


10 


CN 


flux_BGv2 


der.BGv2 


LDC_BGv2 


15000-30000 


1.75-4.75 


1.0 


2 


GS98 


flux_BGvlO 


der.BGvlO 


LDC_BGvlO 


15000-30000 


1.75-3.00 


1.0 


10 


GS98 


fiux_BGvlOCN 


der.BGvlOCN 


LDC_BGvlOCN 


f5000-30000 


1.75-3.00 


1.0 


10 


CN 


flux_BLv2 


der.BLv2 


LDC_BLv2 


15000-30000 


1.75-4.75 


0.5 


2 


GS98 


flux_BLvlO 


der.BLvlO 


LDC_BLvfO 


15000-30000 


1.75-3.00 


0.5 


10 


GS98 


flux_BLvlOCN 


der.BLvlOCN 


LDC_BLvlOCN 


15000-30000 


1.75-3.00 


0.5 


10 


CN 


flux_BSv2 


der.BSv2 


LDC_BSv2 


15000-30000 


1.75-4.75 


0.2 


2 


GS98 


flux_BSvlO 


der.BSvlO 


LDC_BSvlO 


15000-30000 


1.75-3.00 


0.2 


10 


GS98 


ffux_BSvfOCN 


der.BSvlOCN 


LDC_BSvlOCN 


15000-30000 


1.75-3.00 


0.2 


10 


CN 


flux_BTv2 


der.BTv2 


LDC_BTv2 


15000-30000 


1.75-4.75 


0.1 


2 


GS98 


ffux_BTvlO 


der.BTvlO 


LDC_BTvlO 


15000-30000 


1.75-3.00 


0.1 


10 


GS98 


flux_BTvlOCN 


der.BTvlOCN 


LDC_BTvfOCN 


15000-30000 


1.75-3.00 


0.1 


10 


CN 


flux_BZv2 


der.BZv2 


LDC_BZv2 


15000-30000 


1.75-4.75 


0.0 


2 


GS98 


ffux_BZvlO 


der.BZvlO 


LDC_BZvlO 


15000-30000 


1.75-3.00 


0.0 


10 


GS98 



Notes. GS98 - chemical mixture by Grevesse & Sauval ( il998b . CN - chemical mixture enriched in helium and nitrogen, and depleted in carbon. 
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Fig. 1. The NLTE flux derivatives over log Teff (the left panel) and log^ (the right panel) as a function of temperature for the 
Stromgren uy passbands. Three values of \ogg = 3.5,4.0,4.5 were considered. The values of the atmospheric metaUicity and 
microturbulent velocity were assumed as (Z/Zo)atm=l and ^,-2 km/s, respectively. 



scribed above. Subsequently, for each disc position, fi, we inte- 
grated intensities over the uvby and UBVRIJHK bands 



/ IiA,ix)S(A)dA 
/ S(A)dA 



(9) 



where Ix(j-i) is the specific intensity in the passband x and I(A, fx) 
is the monochromatic specific intensity obtained by SYNSPEC. 
The response function, S(A), was interpolated to the wave- 
lengths of computed monochromatic specific intensity by the cu- 
bic spline method. 



The non-linear limb darkening coefficients, a^, were deter- 
mined using the least squares method by minimizing 



20 



X^(aua2, 03, 04) = ^ (.hoCi - hf , 



(10) 



where IldCi is the specific intensity computed with the Claret 
nonlinear law and /, is the corresponding model intensity, at the 
point /I/. 

Similarly to Heyrovsky (120071) . we evaluated the quality of 
our fits by the relative residuals 



20 

r2 = ii 

N 



(11) 
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Fig. 2. A comparison of NLTE and LTE flux derivatives as a function of the central wavelengths of the U BVRIJHK passbands, for 
a models with Teff =20000 K, \ogg - 4.0, {ZIZq)^h^-\, ^f=2 km/s. The left and right panel correspond to the temperature (a-j) and 
gravity (o-J) derivatives, respectively. 
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Fig. 3. Values of Qj (the left panel) and a| (the right panel) as a function of temperature for LTE and NLTE models and the 
Stromgren uy passbands. The remaining parameters of the model are: \ogg - 4.0, {ZIZq).^h^-\, ^,-2 km/s. 



whereas, the conservation of the flux was controlled by comput- 
ing the relative flux excess 



1 I 

r I_,fi dfl- f h,LDClJ^ d/i 

AT, ^ 



(12) 



The limb darkening law, hijj), given in Eq.|6] is defined as 
(Daszynska-Daszkiewicz et al. 2002) 



hx(p) — 2^ _ = 27r 



4 

k=l 



(13) 



1 - y —a" 



Names of the LDC files are given Table[T|and they are coded 
in the same way as the flux and flux derivatives tables (cf. 



Sect. 3.1). Each file contains the following columns: names of 
LDC, Ok, effective temperature, Teff, gravity, log^, metallicity, 
(Z/Z0)atm, mictroturbulent velocity, and values of LDC in the 
Stromgren and Johnson-Cousins-Glass photometric passbands, 
in the order: uvbyU BVRIJHK. 

In Fig. 4 we show distributions of the normalized specific 
intensity in the UBVRIJHK passbands for a model with the 
following parameters: Teff =20000 K, log g=4.0, ^,=2 km/s, 
(Z/ZQ)atm=l- In this figure, we compare the specific intensity 
computed by means of SYNSPEC for 20 equally separated 
points of yU (squares) and the fitted limb darkening function de- 
fined by Eq.|8](solid lines). Intensities for passbands BVRIJHK 
were shifted downward by n ■ 0.1, where n=l,2,...,7. As we can 
see a quality of the limb darkening fit is very accurate. 

To compare values of the NLTE and LTE intensities, in Fig. 5 
we plot Ix [erg-cm"^-s"'-ster"'] as a function of the angle ji for 
the Stromgren uvby passbands. The same model as in Fig. 4 was 
considered. The NLTE intensities get lower values, in particular. 
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Fig. 4. The angular distribution of the normalized specific inten- 
sity for the UBVRIJHK passbands and the NLTE model with 
refF=20000 K, logg=4.0, ^,=2 km/s and (Z/Zq)^^^!- The ac- 
tual model intensities are marked as squares and intensities com- 
puted from LDC as solid lines. For clarity, passbands from B are 
shifted downward by 0.1. 
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Fig. 5. The intensity distribution for the uvby passbands for 
the NLTE (solid lines) and LTE (dashed lines) models with 
Jeff =20000 K, logg=4.0, f,=2 km/s and (Z/Zo)atm=l. 



for the u passband. This is caused by the location of the u filter 
on the B aimer continuum where the difference in the amount of 
energy radiated in LTE and NLTE models is more pronounced 
than in a region to the right from the Balmer jump, where other 
filters (vby) are defined. 

In Fig. 6 and 7, we plot the NLTE limb darkening coeffi- 
cients, a[, 02,02, 04, as a function of Teff in the Stromgren u 
passband, for different values of gravity, \ogg, and metallic- 
ity, (Z/Z0)atm, respectively. All panels have the same scale. As 
we can see, there is a strong dependence of LDC on effective 
temperature, gravity and metallicity. The sensitivity on T^ff is 
stronger for lower values of log ^ = 3.5. The sensitivity to metal- 
licity is similar for all LDC. The effect of the microturbulent 
velocity, ^t, not shown here, is comparable to the effect of metal- 
licity. 

Let us now discuss the accuracy of our determinations of the 
Claret limb-darkening coefficients for the BSTAR2006 NLTE 
models. In Table. 2, we give the average and maximum values 



Table 2. Comparison of our limb-darkening fit quality with 
results of Heyrovsky (2007) for the Claret non-linear law. In 
columns 2 and 3, we give the average and maximum values of 
the relative rms residual, cr, respectively. Columns 4 and 5 con- 
tain the average and maximum values of the relative flux excess, 
|^|. Passbands for which these goodness-of-fit quantities were 
evaluated are given in notes. 



Method 


Average cr 


Max. cr 


Average |^| 


Max. 




% 


% 






20-point° 


0.1 f 9 


0.432 


14.5 X 10"= 


9.10 X IQ-'* 


20-point* 


0.089 


0.251 


9.50 X 10-"^ 


5.28 X 10-"* 


17-point'' 


0.190 


0.625 


8.27 X 10-^ 


5.05 X lO-'* 


1 f -point'' 


0.175 


0.558 


4.57 X 10-5 


2.17 X 10-* 



Notes. " values evaluated for all 12 passbands: uvbyUBVRIJHK. * val- 
ues evaluated for the BV/?/ passbands. ^ values from Heyrovsky J2007t 
for the BVRI passbands. 



of the merit functions, i.e. the relative residuals, cr, and relative 
flux excess, |^|. These values were calculated for the whole 
range of T^s, i.e. (15000,30000) K, and logg, i.e. (1.75, 4.75) 
dex, assuming the solar metallicity, (Z/Z0)atm=l, and the mi- 
croturbulent velocity of ^,-2 km/s. We used 20 equally spaced 
points of yU instead of 17 unequally spaced points as used in LTE 
Kurucz models. Moreover, our NLTE values of I{fj) were com- 
puted closer to the stellar limb, i.e. up to - 0.001, whereas 
the lowest angle of yU in LTE models is 0.01. We considered two 
sets of passbands: uvbyUBVRIJHK and BVRI. The second set 
was used to compare our results with those of Heyrovsky (2007]) 
who made computations for LTE Kurucz atmosphere models us- 
ing 17 and 1 1 points of ju. As we can see the quality of our fit is 
very good. The average value of cr amounts to 0.119% and for 
the worst fitted profile we got cr-QA32%. Our average value of 
cr is smaller than the Heyrovsky's one by a factor of two. When 
we limited our analysis to the BVRI passbans, as in Heyrovsky 
(120071) . the result is even better. In this case, the average and 
maximum values of cr amount to 0.089% and 0.251%, respec- 
tively. This indicates that with the larger number of equally sep- 
arated points of fi in the fitting procedure one reproduces more 
accurately the model intensities. 

Finally, let us check the conservation of the flux. The average 
value of I ^ I from our fitting for NLTE models is slightly worse 
than in Heyrovsky (2007) for LTE models atmospheres, but of 
the same order of magnitude. The maximum values of |^| are 
9.10 X 10 for all 12 passbnds and 5.28 x 10 for BVRI. The 
corresponding average values are 14.5 x 10"^ and 9.50 x 10"^, 
respectively. For larger values of gravities, the fluxes computed 
with the Claret limb darkening law are overestimated whereas 
for smaller values of log g, fluxes are underestimated. This is 
caused by a steeper slope of the intensity near the limb for lower 
values of log g. Consequently, the fitted limb darkening for lower 
gravities is slightly below model intensities. 

4. Uncertainties in photometric pulsational 
observables from model atmospheres 

The atmospheric metallicity, (Z/Zoj^i^, microturbulent velocity, 
and effects of NLTE are the most important factors which can 
affect the photometric amplitudes and phases of a pulsating star. 

In all comparisons, we used a reference model computed 
with the mass of M = IOM0, the OPAL opacity tables, the 
AGSS09 mixture, hydrogen abundance of X-0.1, metal abun- 
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Fig. 6. The NLTE limb darkening coefficients in the u passband as a function of Teff for three values of gravity: log = 3.5 (solid 
line), log = 4.0 (dashed line) and log - 4.5 (dotted line). The solar metallicity, (Z/Zo)atin=l, and microturbulent velocity of ^,=2 
km/s were assumed. 



dance of Z-Q.Q2, without overshooting from a convective core, 
aov=0.0, and NLTE-TLUSTY models of stellar atmospheres 
with the metallicity of (Z/Z0)atin=l-O and the microturbulent ve- 
locity of ^,-2 km/s. 

In Fig. [8] we compare the photometric observables computed 
with LTE and NLTE model atmospheres at the same values of 
metallicity, (Z/Z0)atm - 1, and microturbulent velocity, - 2 
km/s, for the lOM© model in the course of its main sequence 
evolution. We considered the Stromgren uy passbands and the 
first three radial modes: pi, p2, P3. In the left panel, we show 
the amplitude ratios, A„/Av, and in the right one the correspond- 
ing phase differences, - (py, as a function of T^ff. As we can 
see, the fundamental mode is most sensitive to the NLTE effects. 
In general, the NLTE values of the amplitude ratio, A„/Ay, are 
smaller than the LTE ones for the hotter models and larger for 
more evolved models. 

Subsequently we checked the NLTE effect for nonradial 
modes. In Fig.|9l we plot the same photometric observables as 
in Fig.[8] but for the t - I mode considering one pressure mode, 
pi, and one high-order gravity mode, gi4. Fig. [TOl illustrates the 
same but for two { - 2 modes: pi and g22. In the case of non- 
radial modes, the amplitude ratios, Au/Ay, have larger values for 
NLTE computations. We can see also a different behavior and 
values of Au/Ay and (/;„ - ipy for pressure and high-order gravity 
modes, especially for the f = 2 modes. 



Effect of the atmospheric metallicity on the photometric ob- 
servables for the radial modes are shown in Fig.[TT| Here, we 
considered (Z/Z0)atm = 0. 1 in addition to our standard value of 
(Z/ZQ)atm = 1-0. As we can see, the effect of the atmospheric 
metallicity is comparable to the effect of the departure from LTE 
in stellar atmospheres. In the case of nonradial modes the effect 
of (Z/Zo)atm is negligible. 

Influence of the microturbulent velocity, is of the same or- 
der as (Z/ZQ)atm, as has been shown by Szewczuk & Daszynska- 
Daszkiewicz (2010). 

5. Uncertainties in pulsational photometric 
observables from pulsation theory 

Linear nonadiabatic theory of stellar pulsation provides eigen- 
frequencies, corresponding eigenfunctions and information on 
mode excitation. The value of the flux eigenfunction at the level 
of the photosphere enters the expression for the light variation 
and it is called the /-parameter, as it has been introduced in 
Sect. 2. For a given mode frequency, this parameter depends, 
in general, on mean stellar parameters, chemical composition, 
opacity data, subphotospheric convection etc. Moreover, in the 
case of high-order gravity modes, the /-parameter depends also 
on the mode degree, {! (e.g. Daszynska-Daszkiewicz & Walczak 
2010). There are also large differences in values of / between 
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Fig. 7. The NLTE limb darkening coefficients in the u passband as a function of Jefr for six values of the atmospheric metalhcity: 
(Z/Zo)atm=2.0 (solid black line), 1.0 (dashed black line), 0.5 (dotted black Une), 0.2 (sohd grey hne), 0.1 (dashed grey hne) and 0.0 
(dotted grey line). The gravity of logg=4.0 and microturbulent velocity of ^i=2 km/s were assumed. 




Fig. 8. A comparison of the NLTE and LTE values of the amplitude ratios, Au/Ay (the left panel) and phase differences, <fiu - Vy 
(the right panel) for 10 M© main sequence models as a function of T^g for the first three radial modes: pi, p2, P3. The atmospheric 
metalhcity of (Z/Zo).^^^ = 1 and microturbulent velocity of = 2 km/s were assumed. Linear nonadiabatic pulsations were 
computed with hydrogen and metal abundance of X = 0.7 and Z = 0.02, respectively, the OPAL opacities and the AGSS09 
chemical mixture. 
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Fig. 9. The same as in Fig.[8]but for two i - I modes; pi and gi 
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Fig. 10. The same as in Fig.|8]for for two i = 2 modes: pi and g22- 
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Fig. 11. Effect of the atmospheric metallicity on the amplitude ratios, A„/Av (the left panel) and phase differences, ^„ - (fy (the right 
panel) for 10 Mq main sequence models as a function of T^ff for three radial modes: pi, p2, ps. NLTE model atmospheres with the 
microturbulent velocity of ^, = 2 km/s and two values of (Z/Zq)^^ were used. 
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Fig. 12. Effects of the internal abundances of hydrogen, X, and metals, Z, on photometric observables for the first three radial modes 
as a function of T^g for the 10 M© main sequence models. The NLTE models with the atmospheric metallicity of (Z/Z0)atm = 1 and 
microturbulent velocity of ^, = 2 km/s were used. Linear nonadiabatic pulsations were computed with OPAL opacity tables and 
AGSS09 chemical mixture. 



high and low frequency modes of a given degree, f. This is why 
the amplitude ratios and phase differences for p- and high-order 
g-modes behave so different, which is exemplified most pro- 
nounced in Fig. 10. 

Pulsations of the B-type stars are driven by the k mecha- 
nism operating in the „Z-bump" layer The nonadiabatic com- 
plex parameter, /, is also defined in this subphotospheric layer, 
therefore its value strongly depends on metal abundances, 
chemical mixture and hence opacities (Daszynska-Daszkiewicz, 
Dziembowski, Pamyatnykh 2005). 

In Fig. [121 we show effects of the internal abundances of 
hydrogen, X, and metals, Z, on the photometric amplitude ra- 
tios, Au/Ay, and phase differences, (fu - (fy, for the first three 
radial modes: pi, p2, ps. The same model as in Sect. 4 was se- 
lected. Here, we compare computations obtained with Z=0.02 
vs. Z=0.03 and X=0.7 vs. X=0.75. As we can see, the amplitude 
ratios computed with Z=0.03 are larger than those obtained with 
Z=0.02, whereas increasing hydrogen abundance, X, decreases 
the amplitude ratios. The value of i/j„ - (fy changes by about 0.1 
rad at most. These effects are particulary distinct for the radial 
fundamental mode, pi. Comparing Fig.fT2lwith Fig. 8 and 11, we 
can see also that the effects of X and Z on the photometric ob- 
servables are much more pronounced than the effects of NLTE 
and (Z/Z0)atm- Please note that for a given mode degree, {, scales 
of the amplitude ratios and phases differences are the same in all 
figures. 

Effects of the hydrogen and metallicity abundance for the 
£ - 1, pi and gi4 modes is shown in Fig. [13] and for the { - 
2, pi and g22 modes in Fig. [14] As one can see, in the case of 
the nonradial modes these parameters have smaller influence on 
the photometric observables than in the case of radial modes. 
The amplitude ratio, Au/Ay, for the { - 1 mode increases with 
increasing abundances of both Z and X. This is true also for the 
{ - 2, pi mode, whereas the amplitude ratio of the ( - 2, g22 
mode does not behave monotonically. The important result is 
that in the case of nonradial modes the effects of X and Z are 
considerable smaller than the NLTE effects (cf. Fig. 9 and 10), 
which is opposite to the radial mode case. 

Finally, we evaluated effects of the opacity data and chem- 
ical mixture. In Fig. [15] we show the photometric observables 



for the radial modes computed with two sources of the opacity 
tables and two chemical mixtures. Here, we compared compu- 
tations obtained with the OPAL vs. OP tables and the GN93 vs. 
AGSS09 mixture. Again, the largest effects are for the pi mode 
in both the amplitude ratios and phase differences. The ampli- 
tude ratios computed with GN93 are smaller than those obtained 
with AGSS09, because of relatively higher abundance of iron in 
AGSS09. Also computations with the OP data give smaller val- 
ues of A„/Ay comparing to those obtained with the OPAL data. 
Moreover, the maximum of the OP amplitude ratio is shifted to- 
wards higher effective temperature. This is connected with the 
location of the Z-bump layer, which occurs at slightly higher 
temperature in the OP data. A usage of different opacity data 
changes the value of ^„ - (fy by about 0.2 rad at most. The chem- 
ical mixture has rather a minor effect on phase differences. We 
can see that also these input data affect the photometric observ- 
ables of radial modes far more than the atmospheric parameters 
(cf. Fig. 8 and 11). 

Effects of the opacity tables and chemical mixture for non- 
radial modes with { - \ and 2 are shown in Fig. [16] and Fig.fTT] 
respectively. These data affect the ( - \ modes very subtle and 
are far less important than the NLTE effects (cf. Fig. 9). In the 
case of the ( - 2 modes effects of opacities and mixture are 
comparable to the NLTE effects (cf. Fig. 10). 

As we could see, generally, in the case of radial modes ef- 
fects of pulsational parameters on photometric observables are 
much larger than effects of atmospheric parameters discussed in 
previous section. The opposite is true for nonradial modes which 
are more sensitive to the NLTE effects, in particular in the case 
of dipole modes. 



6. Conclusions 

We have presented a comprehensive overview of possible 
sources of uncertainties in theoretical values of the photomet- 
ric amplitudes and phases of early B-type pulsators. These data 
are of particular importance because they serve as tools for mode 
identification in the case of main sequence pulsators, as well as 
contain information on stellar physics. 
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Fig. 15. Effects of the opacity tables and chemical mixture on photometric observables for the first three radial modes as a function 
of Teff for the 10 main sequence models. The NLTE models with the metallicity of (Z/Z0)atm = 1 and microturbulent velocity 
of ^, = 2 km/s were used. Linear nonadiabatic pulsations were computed with X - 0.7 and Z - 0.02. 
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Fig. 17. The same as in Fig.[T5]but for two £ - 2 modes: pi and g22- 



The uncertainties are embedded in stellar model atmospheres 
and nonadiabatic theory of stellar pulsation. The atmospheric in- 
put consists of the flux derivatives over Teft and log g, and limb 
darkening and its derivatives, in photometric passbands. These 
quantities are sensitive to the atmospheric metallicity and mi- 
croturbulent velocity, as well as to the departure from the LTE 
approximation. From pulsation computations, we get the nona- 
diabatic complex parameter, /, whose value depends on chem- 
ical composition, opacities and subphotospheric convection if 
present. If the /-parameter can be derived from observations 
then we get an extra seismic probe by means of which we can 
test this input physics. 

We begun with computations of tables with various quanti- 
ties, needed to evaluate the light variation in a given photometric 
band, for NLTE-TLUSTY model atmospheres. These tables in- 
clude data on the passband flux, flux derivatives over effective 
temperature and gravity, and the non-linear limb darkening co- 
efficients in 12 most popular passbands: uvbyUBVRIJHK. All 
these data are public available and can be retrieved from the 
Wroclaw HELAS web page. 

Then, we studied effects of these input parameters on the 
photometric observables, i.e. amplitude ratios and phase differ- 



ences. We considered the 10 main sequence models and the 
low degree modes with £ - 0, 1,2. In particular, the effect of 
NLTE model atmospheres was studied for the first time. In the 
case of radial modes, this effect is comparable to an influence 
of the atmospheric metalicity and microturbulent velocity. The 
photometric observables of nonradial modes appeared more sen- 
sitive to the departure from LTE. Subsequently, we drew a par- 
allel to effects of parameters related to pulsation computations. 
These comparisons showed that the photometric observables of 
the radial modes are by far more sensitive to the pulsation input. 
In turn, in the case of nonradial modes the NLTE effect became 
more important, especially for the € - I modes. 

A complete seismic model should reproduce not only pulsa- 
tional frequencies but also the observed values of the photomet- 
ric amplitude and phases, which can be translated into the em- 
pirical values of the /-parameter. Our studies showed that such 
comprehensive approach should take into account all inaccura- 
cies in the photometric amplitudes and phases. 
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